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In this paper wc consider how non-Gaussianity of the primordial density perturbation and the 
amplitude of gravitational waves from inflation can be used to determine parameters of the curvaton 
scenario for the origin of structure. We show that in the simplest quadratic model, where the 
curvaton evolves as a free scalar field, measurement of the bispectrum relative to the power spectrum, 
/nl, and the tensor-to-scalar ratio can determine both the expectation value of the curvaton field 
during inflation and its dirricnsionless decay rate relative to the curvaton mass. We show how 
these predictions arc altered by the introduction of self-interactions, in models where higher-order 
corrections are determined by a characteristic mass scale and discuss how additional information 
about primordial non-Gaussianity and scale dependence may constrain curvaton interactions. 

I. INTRODUCTION 

Inflation solves the horizon problem, the flatness problem and the monopole problem. Furthermore, it gives a simple 
way to source primordial perturbations from quantum vacuum fluctuations. Any light scalar field during a period of 
inflation with an almost constant Hubble expansion acquires an almost scale-invariant power spectrum of fluctuations 
that could be the origin of primordial density perturbations [1, 2]. 

The curvaton is one such field which is only weakly coupled and hence decays on a time-scale much longer than 
the duration of inflation [3-7]. Its lightness enables the field to acquire super-Hubble perturbations from vacuum 
fluctuations during inflation. When it decays into radiation some time after inflation has ended, its decay can source 
the perturbations in the radiation density of the universe, and all other species in thermal equilibrium [8, 9]. 

One of the distinctive predictions of the curvaton scenario for the origin of structure is the possibility of non- 
Gaussianity in the distribution of the primordial density perturbations [10-12]. Treating the curvaton as a pressureless 
fluid one can estimate the resulting non-Gaussianity either analytically by treating the decay of the curvaton as 
instantaneous [9 11], or numerically [13, 14], showing that the non-Gaussianity parameter /nl becomes large when 
the curvaton density at the decay time becomes small. 

The non-linear evolution of the field before it decays can also contribute to the non-Gaussianity of the final density 
perturbation. The authors of [15-20] look at the effect of polynomial corrections to the quadratic curvaton potential. 
In some cases the curvaton density can be significantly subdominant at decay and still yield small /nl [17]. For small 
values of /nl, the non-Gaussianity can instead be probed by the trispectrum parameter, jinl- 

Primordial gravitational waves on super-Hubble scales are also present since they are an inevitable byproduct at 
some level of an inflationary expansion. Non-Gaussianity alone could distinguish between the curvaton scenario and 
the conventional inflaton scenarios for the origin of structure since a single inflaton field is not capable of sourcing 
significant non-Gaussianity [21]. But non-Gaussianity and gravitational waves together can give tight constraints on 
curvaton model parameters. Nakayama et al [22] recently studied the effects of the entropy released by the decay 
of a curvaton field with a quadratic potential on the spectrum of gravitational waves that are already sub-horizon 
scale at the decay and consider the possibilities of future direct detection experiments, such as BEG or DECIGO, to 
constrain the parameter space. We will restrict our attention to gravitational waves on super-Hubble scales when the 
curvaton decays which are not affected by the decay, and consider self-interactions of the curvaton field in addition to 
the quadratic potential [20, 23]. This includes scales which contribute to the observed CMB anisotropics, where the 
power in gravitational waves is typically given by the tensor-to-scalar ratio for the primordial metric perturbations. 

In this paper we will investigate how non-Gaussianity and gravitational waves provide constraints on curvaton model 
parameters. For any value of the curvaton model parameters we can obtain the observed amplitude of primordial 
density perturbations on large scales by adjusting the Hubble scale of inflation, which we assume to be an independent 
parameter in the curvaton model. However observational constraints on the tensor-to-scalar ratio places an upper 
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bound on the inflationary Hubble scale, while non-Gaussianity constrains the remaining model parameters. 

We numerically solve the evolution of the curvaton field in a homogeneous radiation-dominated era after inflation 
allowing for non-linear evolution of the curvaton field due to both explicit self-interaction terms in the potential 
and the self-gravity of the curvaton. In particular we consider quadratic and non-quadratic potentials which reduce 
to a quadratic potential about the minimum with self-interaction terms governed by a characteristic mass scale, 
corresponding to cosine or hyperbolic-cosine potentials. Cosine potentials arise for PNGB axion fields and are often 
considered as candidate curvaton fields [24-28]. The hyperbolic cosine is representative of a potential where self- 
interaction terms become large beyond a characteristic scale. In each case we show how the non-linearity parameter 
/nl and tensor-to-scalar ratio, vt, can be used to determine model parameters. 

In Section II we review the perturbations generated during inflation and how these are transfered to the primordial 
density perturbation in the curvaton scenario. In Section III we present the numerical results of our study for three 
different curvaton potentials. We conclude in Section IV. 



II. INFLATIONARY PERTURBATIONS IN THE CURVATON SCENARIO 



In the curvaton scenario, initial quantum fluctuations in the curvaton field, x, during a period of inflation at very 
early times give rise to the primordial density perturbation in the subsequent radiation-dominated universe some time 
after inflation and after the curvaton field has decayed into radiation, e.g., the density perturbation in the epoch of 
primordial nucleosynthesis. This primordial density perturbation is conveniently characterised by the gauge-invariant 
variable, Q, corresponding to the curvature perturbation on uniform-density hypersurfaces [29]. 

Throughout this paper we will use the 5N formalism [11, 30-32] to compute the primordial density perturbation 
in terms of the perturbation in the local integrated expansion, TV, from an initial spatially-flat hypersurface during 
inflation, to a uniform-density hypersurface in the radiation-dominated era 

(i = 5N = N'6x. + \N"5xl+ + ... (1) 

where 5x* = X* ~ (x*) and primes denote derivatives with respect to x*; the local value of the curvaton during 

inflation. 

Quantum fluctuations of a weakly-coupled fleld on super-Hubble scales (/c/a <C H) during slow-roll inflation is well 
described by a Gaussian random field with two-point function 

(Xfe.Xfe.) = {'^^fPx{ki)5\k, + fc2) . (2) 
We define the dimensionless power spectrum V^ik) as 

27r2^ 

The power spectrum of curvature perturbations is thus given, at leading order, by 

= N"Vs^{k) . (4) 



^x(fc) = T^Pxik) (3) 



and we define the spectral index as 



_ rflnlnt; - 1| 
dink ■ 

The connected higher-order correlation functions are suppressed for a weakly-coupled scalar field during slow-roll 
inflation, but non-linearities in the dependence of A'' and hence C on the initial curvaton value in Eq. (1) can lead to 
significant non-Gaussianity of the higher-order correlation functions, in particular the bispectrum 

iCkSk^Cks) = i2nfBi{kuk2,ks)S'ih+k2 + ks). (7) 
The bispectrum is commonly expressed in terms of the dimensionless non-linearity parameter, /nl, such that 

Bc{ki,k2,k3) = lhi^[Pc{k,)Pc{k2) + Pdki)Pc{k3) + P^{k2)Pcik3)] (8) 
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and the running of the spectral index as 
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If the initial field perturbations, 5x*, correspond to a Gaussian random field then it follows from Eq. (1) that /nl 
is independent of the wavenumbers, ki, and is given by 

5 N" 

/nl = g-^- (9) 

In practice non-linear evolution of the field can lead to non-Gaussianity of the field perturbations on large scales and 
a weak scale dependence of /nl [33-35] . 

Current bounds from the CMB on local-type non-Gaussianity require — 10 < /nl < 74 [36]. Large-scale structure 
surveys lead to similar bounds [37]. 



A. Isocurvature field perturbations during inflation 



Perturbations of an isocurvature field, whose fluctuations have negligible effect on the total energy density, can be 
evolved in an unperturbed FRW background and obey the wave equation 

a? 



5x + ^H5x+[^+ml]8x = G, (10) 



where the effective mass-squared is given by = d^V/dx^. During any period of accelerated expansion quantum 
vacuum fluctuations on small sub-Hubble scales (comoving wavenumber k > aH) are swept up to super-Hubble scales 
{k < aH). For a light scalar field, X) with effective mass much less than the Hubble rate during inflation (m^^ <C Hi) 
the power spectrum of fluctuations at Hubble exit is given by 



P^. ^ ( ^ ) for fc = a^H, . (11) 



^^x=:^.'Px--'2^ + '2Vx- (14) 



On super-Hubble scales the spatial gradients can be neglected and the overdamped evolution (10) for a light field is 

given by 

H-'dx^-VxSx- (12) 
where we define the dimensionless mass parameter 

(13) 

Combined with the time-dependence of the Hubble rate in Eq. (11), given by the slow-roll parameter e = —H/H^, 
this leads to a scale-dependence at any given time of the field fluctuations on super-Hubble scales [5, 38] 

d 

dink' 

which is small during slow-roll inflation, e ^ 1, for light fields with \rj^\ <C 1. 

Self-interaction terms in the curvaton potential during inflation only modify the predictions for the power spectrum 
and spectral tilt beyond these leading order results in the slow-roll approximation. However they do lead to time- 
dependence of the effective mass of the x field, so that the effective mass appearing in the expression for the spectral 
tilt may differ from that when the curvaton oscillates about the minimum of its potential some time after inflation. 
In particular the effective mass-squared during inflation could be negative, leading to a negative tilt, An^ < 0, even 
if e is very small. 

The time-dependence of both e and r]^ 

H-^ri^ 2er,^-^2^ (15) 

H-H ~ -2e(?7^-2e) (16) 

during slow-roll inflation driven by an inflaton field with dimensionless mass 77^ = Vij,^/'6H^ and = 
{d^V/dx^d^)/9H^, gives rise to a running of the spectral index in Eq. (14) [39] 

rflnfc -'^'(~^' + ^'^ + ^x)-2ex0. (17) 

In the following we shall make the usual assumption that the curvaton has no explicit interaction with the infiaton, 
so that = and the running is second-order in slow-roll parameters and expected to be very small. Note, 
however, that in the curvaton scenario the tensor-to-scalar ratio and spectral tilt do not directly constrain the slow- 
roll parameters e and 77^ as in single-inflaton-field inflation, so they could be relatively large. 
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B. Transfer to curvaton density 



In the curvaton scenario, these super-Hubble fluctuations in a weakly- coupled field whose energy density is negligible 

during inflation generates the observed primordial curvature perturbation, after inflation if the curvaton conies to 
contribute a non-negligible fraction of the total energy density after inflation. 

As the curvaton density becomes non-negligible one must include the backreaction of the field fluctuations on the 
spacetime curvature. However on super-Hubble scales, k ^ aH, where spatial gradients and anisotropic shear become 
negligible we can model the non-linear evolution of the field in terms of locally FRW dynamics [40] . In the following 
we will employ this "separate universe" picture [32] and we have 

H2.?fl(v.4xi). ,18) 

where Xl = X + ^Xj ^l, Vl and V^l denote the field, Hubble rate, potential and potential gradient smoothed on 
some intermediate scale {aH)~^ <C i < , and dots denote derivatives with respect to the local proper time. 

Once the Hubble rate drops below the effective mass scale, the long- wavelength modes of the field, xl, oscillate 
about the minimum of the potential. Any scalar field with finite mass has a potential which can be approximated 
by a quadratic sufficiently close to its minimum, and the effective equation of state, averaged over several oscillation 
times, becomes that of a pressureless fluid 

Px = {\mlxl + \xl)^a-\ (19) 

Thus the energy density of the curvaton grows relative to radiation, oc a~^. The curvaton must eventually 
decay if it is to transfer its inhomogeneous density into a perturbation of the radiation density. We assume a slow, 
perturbative decay of the curvaton at a flxed decay rate, F <C m (though we note that oscillating fields can also 
undergo a non-perturbative decay, or partial decay at earlier times [41, 42]). 

In this work we will numerically solve for the evolution of the curvaton field until it begins oscillating and determine 
its subsequent energy density. In order to follow the subsequent evolution and eventual decay of the curvaton density 
on time scales, ^ F"-'^, much longer than the oscillation time, ^ m"'^, we adopt the results of Ref. [43]. 

Once the curvaton field behaves as a pressureless fiuid, one can show that phase-space trajectory is determined by 
the dimensionless parameter [43, 44] 

fW 

r/"H^o"''\ F 



In practice one can only treat the curvaton field as a pressureless fiuid once it has begun to oscillate about the minimum 
of its potential. Taking the density of the curvaton when it begins to oscillate, p^.osc — in Eq. (20), we can 

estimate p as [5] 

where the subscript "ose" denotes the time for which i^osc = "^x ^^'^ ™pi = {8ttG)~^^'^ ~ 2.43x IQ-'^^GeV is the reduced 
Planck mass. Although the actual time when the curvaton begins oscillating is also not precisely defined this need 
not be a problem as Cl^^yH/T is a constant while the curvaton is sub-dominant at early times, since oc a oc 

^1/2 

and H (X for a pressureless fluid in a radiation dominated era, and we simply require Xosc/6™pi ^ 1- 

However, Eq.(20) only estimates p in terms of the curvaton field value when the curvaton starts oscillating and we 
have assumed it has a quadratic potential at this time. More generally, to allow for self-interactions of the curvaton 
field that could lead to non-linear evolution after inflation and could still be significant when the curvaton begins to 
oscillate we define a transfer function for the field Xosc — 9{x*) [13] such that 

^ - 6m|i V r • ^ ' 

in order to relate the density of curvaton at late times, as it oscillates about the minimum of its potential, to the 
value of the curvaton field during inflation, x* • 
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C. Transfer to primordial perturbation 



The amplitude of the resulting primordial curvature perturbation depends both on the perturbation in the curvaton 
density, ^p-^j p-^, and the energy density in the curvaton field when it decays. To first-order in the perturbations we 
write 



C = i?,r^ =R,Z- (23) 



where < i?^ < 1 is a dimensionless efficiency parameter related to the fraction of the total energy density in the 
curvaton field when it decays into radiation. Using the separate universe picture, we take derivatives of the same 
function g{x*) defined in terms of the homogeneous background fields in Eq. (22) to determine the linear density 
perturbation and higher-order perturbations in terms of the field perturbations during inflation. We thus have the 
transfer function for linear curvaton field perturbations during infiation into the primordial curvature perturbation 



C = i?,i^=i?,'^^. (24) 



where primes denote derivatives with respect to x*- 

Modelling the transfer of energy from the curvaton field to the primordial radiation by a sudden decay at a fixed 
value of Hdecay = T gives the transfer parameter [5, 9] 



4ptotal - Px. 



(25) 

decay 



However this expression is of limited use if we want to predict the primordial curvature perturbation in terms of the 
inflationary value of the curvaton fleld and its perturbations because this expression refers to the curvaton density at 
the decay time. The curvaton density changes with time and the decay time is not precisely deflned since the decay 
happens over a flnite period of time around iJ ~ F. 

More generally, the transfer parameter, i?^ in Eq. (23), is a smooth function of the phase-space parameter p deflned 
in Eq. (20). One can determine as a function of p numerically, which gives the analytic approximation [44] 

/ 924 \ "^'^^ 

i?^(p)^l-(^l + — . (26) 

A distinctive feature of the curvaton scenario is the possibility that the primordial curvature perturbation may have 
a signiflcantly non-Gaussian distribution even if the curvaton fleld itself is well described by a Gaussian distribution. 

This is due primarily to the fact that the energy density of a massive; scalar field when it oscillates about the minimum 
of its potential is a quadratic function of the field. Simply assuming a linear transfer (23) from a quadratic curvaton 
density to radiation yields [9] 



^ 3 



and hence a primordial bispectrum of local form [45] characterised by the dimensionless parameter 

/nl = ^. (28) 

This provides a good estimate of the non-Gaussianity for a quadratic curvaton with Gaussian distribution when 
/nl > 1- 

Incorporating the full non-linear transfer for a quadratic curvaton density while assuming the curvaton field has a 
Gaussian distribution at a sudden decay, yields corrections of order unity [10, 11, 13] 

■^^^^4i,~3~^- ^^^^ 

Numerical studies [13, 14] confirm that this sudden-decay formula for /nl(-Rx) represents an excellent approximation 

to the actual exponential decay, cx e^^'- /a?, where we take in Eq. (29) to be the linear transfer efficiency 
deflned by Eq. (23). In particular we flnd the robust result /nl > — 5/4 for any value of R^. 
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More generally, if wc allow for possible non-linear evolution of the local curvaton field after Hubble-exit through 
the function g{x*) defined in Eq. (22), and allow for possible variation of the transfer parameter with the value 
of the curvaton field (but still take the curvaton fluctuations to be Gaussian at Hubble-exit) then we have [13] 



/nl — 



4R^ 



5r 

„/2 



R^' jg/g') - 2R^ 



(30) 



This expression follows directly from Eq. (9) when we take A''' 

If we adopt the sudden-decay approximation for Rx{p) then Eq. (30) reduces to [11] 



2e. 91 

g 



/i 



NL 



AR-y 



1 + 



9"g 



5R^ 



(31) 



D. Metric perturbations during inflation 

In most studies of the curvaton scenario it is assumed that the amplitude of scalar or metric perturbations generated 
during inflation arc completely negligible. Indeed the original motivation for the study of the curvaton was to show 
that it was possible for fluctuations in a field other than the inflaton to completely dominate the primordial curvature 
perturbation. However gravitational waves describe the free oscillations of the metric tensor, independent (at first 
order) of the matter perturbations, and some amplitude of fluctuations on super-Hubble scales is inevitably generated 
during an accelerated expansion. The resulting power spectrum of tensor metric perturbations is given by 

The power spectrum of primordial gravitational waves if they can be observed by future CMB experiments, such 
as CMBPol [46], would give a direct measurement of the energy scale of inflation and hence the Hubble rate, H^. 
In practice the amplitude of gravitational waves is usually expressed relative to the observed primordial curvature 
perturbation as the tensor-to-scalar ratio 
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J^S.lxlO^f^) =0.14x(— ^) . (33) 



Current observational bounds from CMB anisotropics are partially degenerate with bounds on the spectral index and 
dependent on theoretical priors, but can be used give rr < 0.24 [36]. Bounds from the power spectrum of the B-mode 
polarisation of the CMB are less model dependent and require < 0.72 [47]. 

The tensor perturbations are massless and the scale dependence of the spectrum after Hubble-exit (32) is simply 
due to the time dependence of the Hubble rate: 

riT = -2e . (34) 

Thus the tilt of the gravitational wave spectrum on very large scales today gives a direct measurement of the equation 
of state during inflation, w = — 1 -|- (2e/3). 

If inflation is driven by a light inflaton field, (f, this inflaton field also inevitably acquires a spectrum of fluctuations 
during the accelerated expansion, V,p, = {H/2tt)1. These adiabatic field perturbations [48] correspond to a curvature 
perturbation at Hubble-exit during inflation 

hV ^ 1 



The scale-dependence of the tensor spectrum (34) together with the time-dependence of e during inflation, given in 
Eq. (16), leads to a scale dependence of the curvature perturbation from adiabatic perturbations 

nc,-l = -6e + 2r]^ , (36) 

where the dimensionless inflaton mass parameter is r]^ = m^/3H'^. Note that the primordial curvature perturbation 
due to canonical inflaton fleld perturbations is eflfectively Gaussian with |/nl|* ^ 1 suppressed by slow- roll parameters 
[21]. 
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In the presence of a curvaton field, the adiabatic perturbations during inflation represent only a lower bound on 
the primordial curvature perturbation and one should add the uncorrelated contributions to the primordial curvature 
perturbation from both the curvaton field (24) and the inflaton field (35): 

^< = (^)'^« + i^^'- w 

For example, if the spectral tilt of the primordial curvature perturbation from a very light curvaton field (14) is 
— 1 « —0.03 and primarily due to the time-dependence of the Hubble rate during inflation, n^^- — 1 ~ Ut — — 2e, 
then we have 16e « 16 x 0.015 = 0.24 and hence P^, « AVt- Hence Vq^ < for rr < 0.3. 

In the following we will assume e is large enough that the inflaton contribution to the primordial curvature pertur- 
bation can be neglected even if the primordial tensor perturbations are potentially observable. 



III. NUMERICAL RESULTS 



In our numerical analysis we have used the separate universe equations (18) to evolve the local value of xl for 
long-wavelength perturbations of the curvaton field. This incorporates both the non-linear self-interactions included 
in the potential of the curvaton, V{xl)i and non-linearity of the gravitational coupling through the dependence of 
the Hubble expansion rate on the curvaton field kinetic and potential energy density in the Friedmann equation. 

We do not solve for the curvaton field evolution during infiation or during (p)reheating at the end of infiation since 
this would be model dependent. Instead we start the evolution with a radiation density such that the initial Hubble 
rate is much larger than the effective mass of the curvaton, consistent with our assumption that the initial value of 
the curvaton field is effectively the same as its value at the end of inflation, x* • 

We evolve the curvaton until it begins to oscillate in the minimum of its potential and can accurately be described 
as a pressureless fluid, in order to exploit earlier work which used a fluid model to study the linear [44] and non-linear 
[13] transfer of the curvaton perturbation to radiation and hence the primordial curvature perturbation. Thus we 
evolve the curvaton fleld until oc a~^. Note that this may be sometime after the time when H = m-)^ since the 
curvaton potential may have significant non-quadratic corrections at this time. 

We need to be able to determine the dimensionless parameter p defined in Eq. (20) which determines the transfer 
parameter Rx{p)- To do so we identiiy 

/ TYX 

p = y Y^"^ ' ^^^^ 

where 

pp^ = n^{i-n^)'^/*^, (39) 

is constant for a pressureless fluid, x, plus radiation. It is straightforward to check that Eq. (38) coincides with the 
deflnition of p given in Eq. (20), which is evaluated in the early time limit, 0. The advantage of our variable 

Pfw is that it can evaluated at late times, so long as the curvaton decay is negligible, T <^ H, whereas at early times 
the curvaton fleld m ay never actually evolve like a pressureless fluid and we may not have a well-deflned early time 
limit for n^y/H/T. 

In our numerical code following the curvaton field evolution we use Eq. (18) with the rescaled time variable r = mt, 
implicitly setting F = 0, such that 

x" + 3hx'+^=0, (40) 

For a quadratic potential we have Vy^jrr? = x a-iid Vjrr? = and hence the evolution of x{t) is independent 
of m. We evolve the curvaton field from an initial value Xi = X* when Hf = 100^^^. This is consistent with the 
usual assumption that the curvaton is a late-decaying field with F <C m. We are then able to determine pfw{x*) 
which approaches a constant as the curvaton density approaches that of a pressureless fiuid at late times. We then 
obtain the actual parameter p in Eq. (38) for a finite decay rate, by multiplying by a finite value of ^/m/T. Thus the 
parameter p is a function of x* and m/T, but not m and F separately. 
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We use the previously determined [44] transfer function Rx{p) given by Eq. (26). Prom Eq. (11) and (24) we then 

have 



3p 



2lT ) 



(42) 



Normahsing the amphtude of the primordial power spectrum to match the observed value on CMB scales, Vq ^ 
2.5 X 10~^ [36], then fixes the amplitude of vacuum fluctuations of the curvaton field during inflation and hence the 
scale of inflation 



= 9.4 X 10" 



P 



mpi. 



(43) 



or, equivalently, the tensor-scalar ratio 



rr = 72 



p'R^{p)mpi 



(44) 



The non-linearity parameter, /nl, is given by Eq. (30). Note that for tt wc must determine not only p but also 
its first derivative, p' , with respect to the initial field value, x*- For the non-linearity parameter, /nl, we also need 
the second derivative, p", and to describe higher-order non-Gaussianity we would need higher derivatives. In terms 
of the parameter p, Eq. (30) becomes 



/i 



NL 



2R^ 



PP 

nl2 



rIp 



1 



(45) 



A. Quadratic curvaton 



We show the results in Figure 1 and 2 for a quadratic curvaton potential. In this case we are able to compare 
our numerical result against an exact analytic expression while the curvaton density remains negligible during the 
radiation-dominated era. In this case the curvaton field is given by 

^ 25/4r(3/4) (mt)i/4 • ^ > 

where Ji/i{mt) is the Bessel function of the first kind of order 1/4. This has the asymptotic solution x — 
1.023%* cos{mt - 3'K/8)/{mt)^/^, and substituting this into Eq. (20) gives 

We see from Figure 1 that Eq. (47) gives an excellent approximation to the numerical results for x* ^ ^Pi- 

Contour plots arc given in Figure 2 for the non-linearity parameter, /nl, and the inflation Hubble scale, if*, (and 

hence tensor-scalar ratio, Vt) for a non-self-interacting curvaton with a quadratic potential. 

Given that the analytic result for p(x*) given in Eq. (47) is an excellent approximation, except for x* ~ mp\, we 

deduce that Xosc = gix*) defined by Eq. (21) is a linear function g{x*) — V^X*- Thus the non-linearity parameter 

/nl is given in terms of R^ in Eq. (29). We have two regimes for the transfer function Rx{p) given by Eq. (26). 

For X* ^ (T/my^^m-pi we have p ^ 1 and hence R^ ~ 1, while for x* ^ {T/my^^mpi we have p <^ 1 and hence 

R^ ~ 0.924JJ. Thus we find from Eq. (29) 

r -5/4 for X* > (r/m)V^mpi 

/nl — S „„ rFml, f m/ M/4 • ('^^) 

I ^-^V ™^ ^* (r/m)'/*mpi 
Potentially observable levels of non-Gaussianity (5 < /nl < 100) are found in a band of parameter space 

X*«(l-4)xl0i^Gev(^)''\ (49) 
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FIG. 1: Dimensionless curvaton parameter Pfw, defined in Eq. (39) 
three different potentials: quadratic potential (dotted blue line), cosine potential with / 



as a function of initial curvaton field value, x*i for 
lO^^GeV (upper red dashed line) 



and hyperbolic cosine potential with / = 10 Gev (lower green dot-dashed line). For comparison, the solid black line shows 
X«/3mpi, which provides an excellent approximation for x* ^ t^pi- 



The degeneracy between values of x* a-nd T /m which would be consistent with the same value of /nl is broken by 
a measurement of the scalar to tensor-ratio, ry. Substituting the approximation (47) in Eq. (43). We have 

if, ~ 4.7 X 10-*-^, (50) 

This yields two simple expressions for iJ* according to whether p ^ 1 and hence R^^ ~ 1 or p ^ 1 and hence 
R-^ ~ 0.924p. We thus have 

J 4.7 X lO-^x* for X* > {T/my/'^mpi 

^* - I 1.5 X lO-^yiJlgi for X, « (r/m)i/4™pj ' (^1) 

Even a conservative bound on the tensor-scalar ratio such as r-p < 1 thus places important bounds on the curvaton 
model parameters. Firstly there is the model-independent bound on the inflation Hubble scale, H^, < 2.7 x 10^^ GeV. 
In the case of a quadratic curvaton potential this imposes an upper bound on the value of the curvaton during inflation 

X* < 5.7 X 10^^ GeV, (52) 

which is consistent with x* < "^Pi required to use the analytic approximation (47). We also find an upper bound on 
the dimensionless decay rate 

— < 0.023 f — 1 , (53) 

and in any case F < 10~^to. For example, for a TeV mass curvaton [49] we require F < 1 GeV. More generally, 
if we require the curvaton to decay before primordial nucleosynthesis at a temperature of order 1 MeV, we require 
F > Hbbn and hence m > lO'^-ffeBN- On the other hand if the curvaton decays before decoupling of the lightest 
supersymmetric particle at a temperature of order 10 GeV, we require F > 10~^^ GeV and hence m > 10""'^* GeV. 

Bounds on the curvaton decay rate due to gravitational wave bounds were also studied recently in Ref. [22], who 
also considered the case where that curvaton oscillations begin immediately after inflation has ended at H < m. 

We note that bounds on the tensor-scalar ratio rule out large regions of parameter space that would otherwise give 
rise to large non-Gaussianity. 

A simultaneous measurement of primordial non-Gaussianity, /nlj and primordial gravitational waves, r^, for a 
non-self- interacting curvaton field with quadratic potential would determine both the energy scale of inflation, H^,, 
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FIG. 2: Contour plots showing observational predictions for a curvaton field with quadratic potential as a function of the 
dimensionless decay rate, log]^Q(m/r), and the initial value of the curvaton, logj^g(x*/GeV). Top: Contour lines for the non- 
Gaussianity parameter /nl (in blue). The dotted black lines correspond to Eq. (28). Middle: Contour lines for inflationary 
Hubble scale, logi(,(H*/GeV). The plotted contour hnes correspond to H, = lO", lO", 10^^ GeV. The black dotted hnes 
correspond to the 2 limits of Eq. (51). Bottom: Contour lines for both the non-Gaussianity parameter, /nl, (blue thick solid 
line) and tensor-scalar ratio, tt, (red dotted line). 
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and the expectation value of the curvaton, x*- It would also determine the dimensionless decay rate F/to, but not 
the absolute value of the mass and decay rate separately. More optimistically, if the gravitational amplitude was large 
enough to determine the tensor tilt, ut and hence e, the scale dependence of the scalar spectrum would determine 
the curvaton mass: 

ml = SvxH! ^ ^ (rif - 1 - "t) ^;^^f^m|, (54) 

However once e is known then from Eq. (35) we also know the curvature perturbation due to inflaton perturbations 

during slow-roll inflation: = (rr/lBe)?'^. If e < 0.02, as is commonly assumed, then our assumption that the 
inflaton perturbations are negligible is no longer valid for rx ^ 0.3. In this case we need to consider a mixed inflaton- 
curvaton model. This inflaton-curvaton model has a much richer phenomenology [50-54] and we leave the study of 
the combined non-Gaussianity and gravitational wave bounds in this scenario to future work. Otherwise we must 
assume e is sufficiently large that the inflaton-generated perturbations remain negligible. 

B. Self-interacting curvaton 

We have seen that non-linear field evolution due to gravitational back-reaction of the curvaton field with a quadratic 
potential is limited to large initial values %* ~ 'mpi which are incompatible with bounds on the tensor-scalar ratio 
in the curvaton scenario with a quadratic potential. However significant non-linear field evolution may arise from 
self-interactions of the curvaton field, due to deviations from a purely quadratic potential. Polynomial self-interaction 
terms of the form Vint cx^ x" where n > 4 have been shown to have a large effect on observational predictions in some 
regions of parameter space [17, 19, 20]. 

Rather than choose a monomial correction term, we choose a functional form that leads to significant corrections 
at a specified mass scale. In particular we are motivated by axion type potentials where the curvaton field has a 
natural range, /. Thus we consider a cosine-type potential, with a smaller mass effective mass for x* f and a 
hyperbolic-cosine potential which has a much larger mass for x* ~ /■ In both cases the corrections lead to a finite 
range X* ^ / for the initial curvaton field. 

1. Cosine potential 

We consider an axion-type potential for a weakly-broken {7(l)-symmetry[5, 55] 

V{x) = (l - cos lm\' _ ^!!^ + . . . , (55) 

where = M'^ / f"^ <C and we have an additional model parameter corresponding to the mass scale f ^ M 
which determines the relative importance of self-interaction terms at a given curvaton field value. It also determines 
a natural expectation value for the curvaton field, x* /■ In the following we assume / < mp\. 

In Figm'e 1 we show the numerical solution for ppw as a function of x*, corresponding to p for a fixed value of m/T. 
As expected we see that for X* / we recover the analytic result (47) as the potential is effectively quadratic and 
self-interactions have a negligible effect. For larger values of x*, the potential becomes flatter than the corresponding 
quadratic potential and we see that ppw, and hence p, can become much larger than woiild be obtained for a quadratic 
correction. Note that the potential (55) is periodic and we can identify pfw{x* + ''^f /'^) = PFw{T^f /'^ ~ X*)- 

We show numerical predictions for the non-Gaussianity parameter, /nl, and the tensor-scalar ratio, in Figure 3. 
Non-linear evolution of the field becomes important for x* ^ /• In particular we see that an upper bound on the 
tensor-scalar ratio no longer places an upper bound on the decay rate F/m as we approach the top of the potential, 
i.e., as X* tt/. 

Modest, positive values of the non-linearity parameter, 1 < /nl < 10, become possible even if the curvaton 
dominates the energy density when it decays (p > 1) if x* > 2.5/, but we never find very large values of /nl > 100. 
Because g" > in Eq. (31) wc have /nl > —5/4, as in the case of a quadratic potential, and we never find large 
negative values of /nl for a cosine-type potential. 

2. Hyperbolic- CO sine potential 

Non-linearity of the cosine potential (55) yields a flat potential with small effective mass during inflation for 
X* ~ /. To consider the effect of self-interactions leading to a larger effective mass we consider a hyperbolic cosine 
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FIG. 3: Contour plots showing observational predictions for a curvaton field with a cosine type potential. The three plots 
show, from top to bottom, observational parameters for cosine potentials with / = lO^'^GeV, / = lO^^GeV and / = lO^^GeV, 
respectively, as a function of the dimensionless decay rate, log]^Q(m/r), and the initial value of the curvaton, x*/GeV. Thick 
solid blue contour lines show bispectrum amplitude, /nl, decreasing from left to right. Dotted red contour lines show the 
tensor-scalar ratio rr, also decreasing from left to right. 
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potential which becomes an exponential function of the curvaton field at large field values, as may be expected due 
to supergravity corrections. 

V{x) = (cosh - l) ^ \m\- + + . . . . (56) 

In the following we assume / < mpi. As in the case of the cosine potential, this also yields a natural range for x* ^ f ■ 
In the case of a hyperbolic potential, the field becomes heavy relative to the Hubble scale and evolves rapidly for 
values of x* much larger than /. In particular the requirement that the curvaton have an effective mass less than 
0.1_ff at the start of our numerical solutions imposes the constraint x* < 5/. 

In Figure 1 we show the numerical solution for p as a function of x* ^ fixed value of ra/T. As expected we see 
that for X* / we recover the analytic result (47) when the potential is effectively quadratic. However for the a 
hyperbolic potential we see that due to the steeper potential the effective energy density when the curvaton decays, 
determined by the parameter p, becomes less than the quadratic case for x* ^ f ■ 

We show numerical predictions for the non-Gaussianity, /nl, and the inflation Hubble scale, -ff*, (and hence the 
tensor-scalar ratio) in Figure 4. The non-linear correction g"g/g'^ in Eq. (31) becomes negative for x* ~ / and we 
can obtain large negative values of /nl- 

However we find that the bound on the tensor-scalar plays an important role. Regions of parameter space which 
yield large negative /nl also give large tensor-scalar ratios. In regions where 1 and the curvaton dominates when 
it decays we have ^ 1 and both the tensor-scalar ratio and the non-linearity parameters become functions solely 
of X*- III tfiis regime, we have, from Eqs. (44) and (31) 



2 \5'mpi 
4 Vff 



/nl - 7 ( ^ ) - (58) 



which are both clearly functions of x*- Indeed formally we can eliminate ^(x*) and its derivatives in order to write 

/nl ^ -^/H^^Pi (V^)' • (59) 

Hence the contours of equal values of both tt and /nl become horizontal on the right-hand-side of Figure 4. For 
example, with / = 10^^ GeV a weak bound on the tensor-scalar ratio of Tt < 1 requires /nl > —1000. A stronger 
bound < 0.1 requires /nl > —100. 

Of course {^/rr)' is not an observable parameter, but if we assume that -^/tt is a smooth function of x*// we can 
estimate {y/fr)' ~ {y/^)/f and hence 

mpi — 

/NL~--p\A¥- (60) 

This semi-empirical relation appears to hold for sufficiently small T/m and it would be interesting to see if this is also 
the case for the polynomial correction terms [17, 19]. 

Unlike the case of a cosine-type potential we still have a strict \ipper bound on the decay rate, as in the case of a 
purely quadratic potential. Thus, although there are regions of parameter space for x* ^ f i where the non-Gaussianity 
can be small even if the curvaton is subdominant when it decays, p <C 1, we find that these regions correspond to 
large values for the tensor-scalar ratio and are excluded by bounds on primordial gravitational waves. 



IV. DISCUSSION AND CONCLUSIONS 



In this work we have investigated the numerical evolution of a curvaton field from its overdamped regime after 
inflation until it decays into radiation. We have shown how measurement of both the non-linearity parameter, /nl, 
and the tensor-to-scalar ratio, ry, provide complementary constraints on the model parameters. We did this for 
three different curvaton potentials: the quadratic potential, axion-like cosine potentials and hyperbolic potentials. As 
expected both the cosine and the hyperbolic potentials recover the quadratic regime when x* ^ f- 

For the simplest quadratic potential for the curvaton, bounds on the tensor-to-scalar ratio place an upper bound 
on the dimensionless decay rate, ruling out large regions of parameter space that would yield a large primordial 
non-Gaussianity in the distribution of scalar perturbations. Simultaneous measurement of both the non-linearity 
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FIG. 4: Contour plots showing observational predictions for a curvaton field with a hyperbolic-cosine-type potential. The three 
plots show, from top to bottom, observational parameters for potentials with / — lO^^GeV, / — lO^^GeV and / = lO^^GeV, 
respectively, as a function of the dimensionless decay rate, log]^Q(m/r), and the initial value of the curvaton, x*/GeV. Thick 
solid blue contour lines show bispectrum amplitude, /nl, increasing from top to bottom. Dotted red contour lines show the 
tensor-scalar ratio, rr, decreasing from left to right. 
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parameter, /nl, and the tensor-to-scalar ratio, rr, can determine both the expectation value of the field during 
inflation, x*, and the dimensionless decay rate, F/m. 

In the conventional inflaton scenario for the origin of structure we have three free parameters: the inflation scale 
and two slow-roll parameters, e and r]^. These can be determined by power of the primordial scalar perturbations, 
V(^, the tensor perturbations, Vt, and the spectral index of the scalar spectrum, n^. The spectral index of the tensor 
spectrum, if niciasurable, would give a valuable consistency check [46]. Another important consistency condition for 
canonical, slow-roll inflation is that the primordial density perturbations should be Gaussian and the non-linearity 
parameter, /nl, should be much less that unity [21]. 

In the curvaton scenario with a simple quadratic potential we have 5 free parameters: the inflation scale iJ*, the 
expectation value of the curvaton during inflation x*, the decay rate of the curvaton relative to its mass, T/m, and 
the slow roll parameters e and rj^ = ml/3H^. For a curvaton, we find that H^, x* and T/m are determined by 
the primordial scalar perturbations, P^, the tensor perturbations, Vt, and the non-linearity parameter, /nl, but the 
mass and decay rate of the curvaton are not separately determined. The two slow-roll parameters e and r]^ arc then 
determined by the two spectral indices n^, and tit- 

Another natural observable in the curvaton model is the scale dependence of the non-linearity parameter, defined 
as [33] 

_ dln|/NL| 

In the curvaton scenario this is given by a simple expression [28, 35] 

where we define 773 = rrip^V'" /V. This can be rewritten in terms of observable quantities and 7^3 



?^/nl ='73— (63) 
12v2 ,/NL 

Thus it offers the possibility of testing the curvaton self interactions. Future observations may be able to detect 
I/nl'^/nlI > 5 [56], corresponding to \r]s\^/rT > 17. For the quadratic potential we have the consistency condition 
^/nl = 0. 

Deviations from a quadratic potential introduce at least one further model parameter, /, corresponding to the mass 
scale associated with the non-linear corrections. This leads to a degeneracy in model parameters consistent with the 
five observables , Vt, /nl, and tit, but this can be broken by a measurement oi rify^. 

In the case of a cosinc-typc curvaton potential the self interaction corrections became important near the top of 
the potential, i.e., when x* ^ T^f [28] and the tensor-to-scalar ratio no longer places an upper bound on F/m. As 
for a quadratic curvaton, we still find /nl > — 5/4 and hence any large non-Gaussianity, |/nl| ^ 1, has positive /nl- 
But for X* ~ / we have % ~ — (mpi//^) < 0, and if / is well below the Planck scale there could be strong scale 
dependence. 

In the case of a hyperbolic-type potential /nl can become large and negative, for x* ~ /• However the tensor-to- 
scalar ratio again plays an important role, in this case placing a lower bound on /nl, e.g., /nl > —100 for tt < 0.1 
when / = 10^^ GeV. In this regime we find 773 ~ {mpi/f^) > 0, which can be large, leading to strong scale dependence 
for / mpi, with Uf^^^ < for /nl < 0. 

Running of either the scalar tilt, a^, or the non-linearity, a/j^^ [28], yields additional information about the higher 
derivatives of the potential, and in particular curvaton-inflaton interactions which we have assumed are negligible in 
our analysis. 

Significant non-Gaussianity in the primordial perturbations opens up the possibility to extract information from 
the higher-order correlations in the scalar spectrum, such as the trispectrum [13, 16, 19, 57, 58] 

Tcih, k2, ka, ki) = ^5nl [Pc{k2)Pdk3)Pdki) + 3 perms] + ^/^l [Pdki3)Pdk3)Pdki) + H perms] . (64) 

which are sensitive to higher-order derivatives of the expansion history with respect to the curvaton field value during 
inflation through fi(NL = {25/54:)N"' /N'^. Differentiating Eq. (30) we obtain 



25 

5nl = ^ 



Rx^ 9'^ Rx^ V 9'^ 9' J R^^ V 5'" 9'^ 



(65) 
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which using the sudden-decay approximation can be written as [13, 57] 

gjsL find its scale dependence ngj^,^ [34, 35], thus provide additional observable parameters which then offer consistency 
conditions for generalised curvaton models such as the cosine or hyperbolic potentials. In practice we require more 
accurate numerical simulations than those used in this work to reliably determine the required higher-derivatives with 
respect to the initial field value across the range of model parameters used in this paper and we leave this for future 
work. 
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